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1 A hierarchical turbulence model 


The proposal m for a theory of the corrections to the OK theory (“in¬ 
termittency corrections”) is to take into account that the Kolmogorov scale 
itsef should be regarded as a fluctuating variable. 

The OK theory is implied by the assumption, for n large, of zero average 
work due to interactions between wave components with wave length < 
= t n and components with wave length > n K~ n io (t'o being the 
length scale where the energy is input in the fluid and k a scale factor 
to be determined) together with the assumption of independence of the 
distribution of the components with inverse wave length (“momentum”) in 
the shell [n n 1 kk' 1 ]^ 1 , |5) p.420]. 

It is represented by the equalities 


3 

<A = V P+l)i' 
tn A+1 


v = |u|, v 6 R 3 


( 1 . 1 ) 


interpreted as stating an equality up to fluctuations of the velocity com¬ 
ponents of scale K,~ n £o, i.e. of the part of the velocity field which can be 
represented by the Fourier components in a basis of plane waves localized in 
boxes, labeled by i = 1 ,... n 3n , of size n~ n £o into which the fluid (moving 
in a container of linear size £q) is imagined decomposed (a wavelet represen¬ 
tation) so that (n + 1, i') labels a box contained in the box (n, 1 ). 

The length scales are supposed to be separated by a suitably large scale 
factor k (i.e. i n = n~ n £o = n~ 1 l n -i) so that the fluctuations can be consid¬ 
ered independent, however not so large that more than one scalar quantity 
(namely -) suffices to describe the independent components of the ve¬ 
locity (small enough to avoid that “several different temperatures will be 
present among the systems (n + 1 ,j')” inside the containing box labeled 
(n, j), and the v n +ij distribution “will not be Boltzmannian for a constant 
temperature inside”, n p- 2 ]). 

The distribution of v )) +1 ■ is then simply chosen so that the average of the 
v )) +1 ■ is the value v)Kk if the v ^ +1 ■ on scale n + 1 gives a finer description 
of the field in a box named j contained in the box named i of scale larger 
by one unit. 
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Among the distributions with this property is selected the one which 




( 1 . 2 ) 


m =0 i=l 


with Wo a constant that parameterizes the fixed energy input at large scale: 
the motion will be supposed to have a 0 average total velocity at each point; 


hence Wq 3 can be viewed as an imposed average velocity gradient at the 
largest scale £q. 


The Wi n = Wd 3 is then interpreted as a velocity variation on a box of 
scale 4«T n or n~ n as £q will be taken 1. The index i will be often omitted 
as we shall mostly be concerned about a chain of boxes, one per each scale 
n~ n ,7i = 0,1,..., totally ordered by inclusion ( i.e. the box labeled ( i,n ) 
contains the box labeled [i\ n + 1)). 

The distribution of the energy dissipation W n ^ d = v[] i in the hierarchi¬ 
cally arranged sequence of cells is therefore close in spirit to the hierarchical 
models that have been source of ideas and so much impact, at the birth of the 
renormalization group approach to multiscale phenomena, in quantum field 
theory, critical point statistical mechanics, low temperature physics, Fourier 
series convergence to name a few, and to their nonperturbative analysis, 
either phenomenological or mathematically rigorous, 

The present turbulent fluctuations model can therefore be called hier¬ 
archical model for turbulence in the inertial scales. It will be supposed to 
describe the velocity fluctuations at scales n at which the Reynolds number 
is larger than 1, i.e. as long as V " K io > 1. 

The description will of course be approximate, m Sec.3]: for instance 
the correlations of the velocity gradient components are not considered (and 
skewness will still rely on the classic OK theory, [7J Sec.34]). 

Given the distribution (and the initial parameter Wo) it “only” remains 
to study its properties assuming the distribution valid for velocity profiles 
such that v n n~ n £o > v after fixing the value of k in order to match data in 
the literature (as explained in [9l Eq.(12)]). As a first remark the scaling 
corrections proposed in [12] can be rederived. 

Tf the box A = (n,j) C A' = (n — 1 ,j') then the distribution H(W\W^) of 
Wa = is conditioned to be such that ( W) = ac _1 Wa'; therefore the maximum 
entropy condition is that — f Il(W\W') log IV(W\W')dW — Aa / WTUW\W')dW, 
where Aa is a Lagrange multiplier, is maximal under the constraint that ( W) = 
W'k - 1 : this gives the expression, called Boltzmannian in [8], for W(W\W'). 
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The average energy dissipation in a box of scale n can be defined as the 

def 

average of e n = W n £~ n , i n = £qK~ n \ the latter average and its p”th order 
moments can be readily computed to be, for p > 0: 




log ( gn ) 

- login 

, t Wn s P v nTp 


logT(l + p) 
log K 




(1.3) 


Cp — g + t p 


The Wn being interpreted as a velocity variation on a box of scale £q k n , 
the last formula can also be read as exressing the (()p ) ~ r^ p with 
1 _ io g r(g+i) 

^P 3 log AC 

The T p is the intermittency correction to the value I: the latter is the 
standard value of the OK theory in which there is no fluctuation of the 
dissipation per unit time and volume this gives us one free parame- 

<-n 

ter, namely k, to fit experimental data: its value, universal within Ruelle’s 
theory, turns out to be quite large, k ~ 22.75, [3; fitting quite well all 
experimental p- values (p < 18). 

Other universal predictions are possible. In [9j a quantity has been 

studied for which accurate simulations are available. 

If W is a sample (Wq. W\, ...) of the dissipations at scales 0,1,... for the 

distribution in the hierarchical turbulence model, the smallest scale n( W) 
1 

at which WnloK~ n — v occurs is the scale at which the Kolmogorv scale is 

1 

attained (i.e. the Reynolds number Wn J n becomes < 1). 

Taking £q = l,z; = 1, at such (random) Kolmogorov scale the actual 
dissipation is £ = IK n (w) K ”^ W ^ with a probability distribution with density 
P*(0- If w k = \v^ k ± then IT n = Wqw\ ■ ■ ■ vj n and the computation of P*(£) 
can be seen as a problem on extreme events about the value of a product of 
random variables. Hence is is natural that the analysis of P* involves the 
Gumbel distribution <j>(t) (which appears with parameter 3), [9]. 

The P* is a distribution (universal once the value of k has been fixed to 
fit the mentioned intermittency data) which is interesting because it can be 
related to a quantity studied in simulations. 

It has been remarked, [9], that, assuming a symmetric distribution of 

1 

the velocity increments on scale K~ n whose modulus is i , the hierar¬ 
chical turbulence model can be applied to study the distribution of the 
velocity increments: for small velocity increments the calculation can be 
performed very explictly and quantitatively precise results are derived, that 
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can be conceivably checked at least in simulations. The data analysis and 
the (straightfoward) numerical evaluation of the distribution P* is described 


below, following 

2 Data settings 

Let £q,i> = 1 and let W = (Wo, Wi, ...) be a sample chosen with the distri¬ 
bution 



( 2 . 1 ) 


with Wq, k given parameters; and let v = (vq,vi,...) = (W 0 3 , W x 3 ,...). 


Define n(W) = n as the smallest value of i such that W? n 1 = ViK 1 < 
1: n( W) will be called the “dissipation scale” of W. 


Imagine to have a large number J\f of p-distributed samples of W’s. 
Given h > 0 let 


pm = Wwithn(W) = n) n (£ < (W n /W 0 )^ n < £ + hj) (2.2) 


hence hP* (£) is the probability that the dissipation scale n is reached with 

£ in [£,£ + /i]. Then P*(£) d = is the probability density that, at 

the dissipation scale, the velocity gradient is between £ and £ + h. 

The velocity component in a direction is v n cos'd: so that the probability 
that it is in d£ with gradient and that this happens at dissipation scale 
= n is d£ times 



(2.3) 


Let 



(2.4) 


that is the probability distribution of the (normalized radial velocity gradi¬ 


ent) and 



<?m = / h£P(£)£ m (2.5) 
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its momenta. To compare this distribution to experimental data [B] it is 
convenient to define 



We have used the following computational algorithnm to P(£): 

• (1) Build a sample (■i) W« = (W 0 , Wi, ..., W n ,.. .) 

• (2) Stop when n = n* such that > 1 > Wn^K~ n 

• (3) Evaluate raj = int(£i/h ) + 1 where = K ni (W^/Wq) 1 ^ 

• (4) goto to (1) during N times 
Then, the distribution P(C) is given by 



P(mh — h/2) = h l P(m) 


(2.7) 


where x(^4) = 1 if A is true and 0 otherwise. It is convenient to define the 
probability to get a given m value as 



( 2 . 8 ) 


where 5(n,m ) is the Kronecker delta. Once obtained Q(m), we can get 
recursively P(m): 




P(m + 1) = P(m) 


m 


and the momenta distribution is then given by: 



( 2 . 10 ) 


Finally, the error bars of a probability distribution (for instance P) are com¬ 
puted by considering that the probability that in N elements of a sequence 
there are n in the box m is given by the binomial distribution: 



( 2 . 11 ) 
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From it we find 

(' n)m = Np(m ) , ((n — (n) m ) 2 ) m = P(m) (l — P(m)) /N (2.12) 

where (.) = ^2 n .D m (n, N). Therefore, the error estimation for the P prob¬ 
ability is given by 

P(m) ± 3 (P(m) (l - P{rn)) /N) 1/2 (2.13) 

We have done 15 simulations with k = 22, N = 10 12 realizations (10 4 cycles 
of size 10 7 ) and different values of Wq\ 10', 10 s , 5 x 10 8 , 10 9 , 2 x 10 9 , 3 x 10 9 , 
4 x 10 9 , 5 X 10 9 , 10 10 , 2 x 10 10 , 5 x 10 10 , 7 x 10 10 , 10 11 , 2 x 10 11 and 5 x 10 11 . 

1/3 

We also use the Reynold’s number R = W 0 . 


1 
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Figure 1: Distribution of events that reach the Kolmogorov scale n~ n for 
different values of the Reynold’s numbers R and k = 22. The total number 
of events is 10 12 . 

The size of k has been chosen to fit the data for the intermittency expo¬ 
nents C p and it is quite large (k = 22, 0): this has the consequence that the 
Kolmogorov scale is reached at a scale K~ n with n = 2,3 and very seldom 
for higher scales, at the considered Reynolds numbers. That can be seen 
in Figure 1 where we show the obtained distributions of hj, i = 1,..., IV. 
In Figure 2 we see the average value of n and its second momenta. We see 
that for low Reynold’s numbers the values is almost constant equal to 2 and 
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Figure 2: Momenta of the n-distribution. Left: average value. Right: second 
momenta 
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Figure 3: Plot as a function of £ = mh of the logarithm of the probability, 
login Q(£)> with Q(£) = that the Kolmogorov scale is reached at scale 

m = f, for different Reynold’s numbers and h = 10~ 3 (£). 


from R ~ 2000 it begins to grow. The second momenta shows a minimum 
for R ~ 1000 where almost all events are in hi = 2. 
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Figure 4: logiop(z) distribution for different Reynold’s numbers. Central 
figure: a(1259.92) = -0.08, a(1442.25) = -0.12, a(1587.40) = -0.24 and 
a(1709.98) = -0.52. Right figure: a(3684.03) = -0.51, a(4121.29) = -0.55, 
o(4621.59) = -0.58, a(5848.04) = -0.61 and a(7937.01) = -0.64. 


The measured distribution of Q(£) (see Figure 3) reflects the superpo¬ 
sition of two distributions: the values of £ associated to the hi = 2 and 
to hi = 3 events. Moreover, for small values of R the overall distribution 
is dominated by the events hi = 2 and for large Reynold’s numbers it is 
dominated by the hi = 3 events. At each case the form of the distribution 
is different: for small R, log\oQ(C) is quadratic in £ and for large R is linear 
in £. 

The behavior of Q defines the behavior of p(z). In Figure 4 we see the 
p(z) behavior. We again see clearly how for low R values the distribution is 
non sensitive to the values of R and it is Gaussian. For intermediate values 
of R the exponential of a quadratic function is a good fit for the measured 
distribution and large enough values of z but its parameters parameters 
depend on R. Finally for R large of 3000 the distribution changes and its 
behavior for large z values seems to be fitted very well by a linear funcion 
with a A-depending slope. 

It is interesting to show the dependence of the momenta of P. a n , as 
a function of R. In Figure 5 we see their behavior. We can naturally 
identify three regions: Region I (A £ [0,1000]) where the momenta are 
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Figure 5: Measured momenta of the distribution P(£), a n , vs Reynold’s 
number. 


almost constant, Region II (R E [1000,4000]) where the moments grow with 
R and Region III (R E [4000, 00 ]) where relative moments tend to some 
asymptotic value. 

Experimental data can be found in [6j and are illustrated by the two 
plots in Figure 6 taken from the cited work 

which give the function log 10 p(z). i.e. the probability density for observing 
a normalized radial gradient z as a function of z = £/ y/ (£ 2 ) in the case of 
homogeneous isotropic turbulence (HIT) (i.e. Navier-Stokes in a cube with 
periodic boundaries) or in the case of Raleigh-Benard convection (RBC) 
(NS+heat transport in a cylinder with hot bottom and cold top). The 
results of Fig. 6, for z > 0 should be compared with those of Fig.4 at the 
corresponding Reynolds numbers. In both cases we see that the distribution 
for high Reynold numbers have linear-like behavior for large z-values. In fact 
for the HIT case and R = 2243 we can fit a line with slope —0.77 in the 
interval z E [3.3,5.57]. Also in the RBC case we can do a linear fit with 
slope —0.42 {z E 5.14,9.69]) for R = 4648. The value obtained is similar to 
the ones we computed on Figure 4. 

In figure 7, we can compare the measured flatness in our numerical exper- 
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Z=r) V /(3 V ) 
x x ' x x'rms 



z=d v /(3 v ) 

x x ' x x'rms 


Figure 6: Measured p(z) by Schumaher et al. [6] for different Reynold’s 
numbers 


iment with the observed by Schumaher et al. [6]. We see that the values are 
similar for small and large Reynold’s numbers but there is a peacked struc¬ 
ture for intermediate values due to the relevant discontinuity when passing 
from hj = 2 to hi = 3 events. 

All these results shows that the important aspect of the experiments is 
quite well captured with the only paramenter k available for the fits, i.e. a 
strong deviation from Gaussian behavior and the agreement of the location 
of the abscissae of the minima of the tails in the second case at the maximal 
Ho; this feature fails in the first case (HIT) as the abscissa is about 30: is 
it due to a too small Reynolds number?. This seems certainly a factor to 
take into account as the curve appears to become independent of R, hence 
universal as it should on the basis of the theory, for R > 4000. 

In conclusion the results are compatible with the OK theory but show 
important deviations for large fluctuations because the Gumbel distribution 
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Figure 7: Measured flatness ( 04 / 02 ) compared with the results by Schuma- 
her et al. [6j for different Reynold’s numbers 


does not show a Gaussian tail. 

All this has a strong conceptual connotation: the basic idea (ie the pro¬ 
posed hierarchical and scaling distribution of the kinetic energy dissipation 
per unit time) is fundamental. 

The need to assign a value to the scaling parameter n is quite interesting: 
in the renormalization group studies the actual value of k is usually not 
important as long as it is k > 1. Here the value of k is shown to be relevant 
(basically it appears explicitly in the end results and its value ~ 20 must, 
in principle, be fixed by comparison with simulations on fluid turbulence). 

Acknowledgements: The above comments are based on numerical 
calculations first done by P. Garrido and confirmed by G. Gallavotti. This 
is the text of our comments (requested by the organizers) to the talk by 
D. Ruelle at the CHAOS15 conference, Institut Henri Poincare, Paris, May 
26-29, 2015. 
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